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Abstract 

A  fuel  cell  is  a  promising  energy  conversion  system  that  will  eventually  become  the  first-choice  for  producing  power  because  of  its  clean 
or  zero-emission  nature.  A  steady-state,  two-dimensional  mathematical  model  with  pressure  and  phase  change  effects  for  a  single  PEM  fuel 
cell  was  developed  to  illustrate  the  inlet  humidification  and  pressure  effects  on  proton  exchange  membrane  (PEM)  fuel  cell  performance.  This 
model  considers  the  transport  of  species  along  the  channel  as  well  as  water  transfer  through  the  membrane.  It  can  be  used  to  predict  trends  of 
the  following  parameters  along  the  fuel  cell  channels:  mole  number  of  liquid  water  and  water  vapor,  pressure,  temperature,  density,  viscosity, 
velocity,  saturation  pressure,  pressure  drop,  vapor  mole  fraction,  volume  flow  rate,  required  pumping  power  and  current  density. 
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1.  Introduction 

Proton  exchange  membrane  (PEM)  fuel  cells,  using 
hydrogen  as  fuel,  emitting  water  and  operating  at  low  temper¬ 
ature  for  quick  startup,  have  been  widely  recognized  as  the 
most  promising  candidates  for  replacing  the  internal  com¬ 
bustion  engine  in  automobiles,  and  for  replacing  batteries 
in  portable  and  micro  applications.  In  recent  years,  research 
and  development  on  fuel  cells  have  accelerated  and  the  PEM 
fuel  cell  technology  has  been  successfully  demonstrated.  But 
many  key  challenges  affecting  PEM  fuel  cell  performance 
still  exist,  and  much  more  efforts  on  improving  the  perfor¬ 
mance  of  PEM  fuel  cell  become  even  more  crucial  than  ever 
before.  Useful  mathematical  models  can  provide  powerful 
tools  for  the  analysis  and  optimization  of  fuel  cell  perfor¬ 
mance. 

Costamagna  and  Srinivasan  [1]  gave  a  very  good  review 
regarding  fuel  cell  science  and  technology  up  to  the  year 
2000.  Another  recent  review  made  by  Yao  et  al.  [2]  pre¬ 
sented  both  empirical  performance  models  and  theoretical 
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models.  In  the  early  1990s,  the  pioneering  work  on  PEM  fuel 
cell  model  development  was  done  by  Bernardi  and  Verbrugge 
[3,4],  and  Springer  et  al.  [5]  who  formulated  one-dimensional 
and  isothermal  models  for  the  gas-diffusion  layer,  active 
catalyst  layer  and  ion-exchange  membrane.  Only  the  direc¬ 
tion  perpendicular  to  the  membrane  was  considered.  The 
model  employed  water  diffusion  coefficient,  electro-osmotic 
drag  coefficients  and  membrane  conductivities  to  predict  the 
change  of  membrane  resistance  with  current  density.  The 
temperature  was  assumed  to  be  constant  and  these  models 
were  unable  to  simulate  the  flow  behavior  along  the  chan¬ 
nels. 

Compared  with  one-dimensional  model,  a  two- 
dimensional  mathematical  model  is  preferred  for  water 
and  heat  management  analysis,  as  the  temperature,  pressure 
and  water  varies  along  the  channel  as  well  as  across  the 
membrane.  Fuller  and  Newman  [6]  developed  a  non- 
isothermal  model  by  including  material  balances  in  the 
channel,  concentration  and  temperature  gradients  along  the 
channel  as  well  as  across  the  membrane  surface.  In  the 
model  developed  by  Nguyen  and  White  [7],  the  focus  was 
on  the  transport  and  reaction  in  the  MEA  and  heat  balances 
in  the  gas  channels.  Subsequently,  an  improved  model  was 
developed  by  Yi  and  Nguyen  [8]  to  compare  different  fuel 
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Nomenclature 

a 

water  vapor  activity  in  stream 

Across 

the  cross-section  area  of  channel  (cm2) 

A  stack 

the  heat  transfer  area  in  a  control  volume  (cm2) 

CMW 

concentration  of  water  at  interface  of  the  mem¬ 
brane  (mol  cm-3) 

Cp,i 

heat  capacity  of  species  i  (J  mol-1  k-1) 

d 

channel  height  (cm) 

D 

hydraulic  diameter  of  channel  (cm) 

D° 

a  parameter  used  in  the  expression  for  diffusion 
coefficient  of  water  (cm2  s-1) 

Dmw 

diffusion  coefficient  of  water  (cm2  s-1) 

f(x) 

friction  factor 

F 

Faraday  constant,  96487  C  equivalent-1 

h 

channel  width  (cm) 

I 

current  (A) 

Kx) 

current  density  (A  cm-2) 

1° 

exchange  current  density  for  the  oxygen  reac¬ 
tion  (A  cm-2) 

kc 

condensation  rate  constant  (s-1) 

kp 

hydraulic  permeability  of  water  in  the  mem¬ 
brane  (cm2) 

L 

length  of  channel  (cm) 

Mi 

molecular  weight  of  species  i  (g  mol-1) 

equivalent  weight  of  a  dry  membrane 
(g  mol-1) 

nd 

electro-osmotic  drag  coefficient 

nQ 

the  mole  number  of  electrons  needed  per  sec¬ 
ond  for  1  A  of  current  (mol  s-1  A-1) 

N 

mole  number  of  species  in  the  stream  (mol  s  - 1 ) 

NE 

number  of  electrons  (A-1  s-1) 

Nch 

number  of  channel  (s) 

P 

local  pressure  (Pa) 

Pi 

partial  pressure  of  species  i  (Pa) 

dp 

pressure  drop  (Pa) 

P  pump 

pumping  power  (W) 

Q 

volume  flowrate  (m3  s-1) 

Re 

Reynolds  number 

Ru 

universal  ideal  gas  constant 
(8.3144  J  mol-1  K-1) 

RH 

relative  humidity 

T 

temperature  of  stream  (K) 

Ts 

temperature  of  stack  (K) 

u 

overall  heat-transfer  coefficient 
(J  s-1  cm-2  °C_1) 

V 

flow  velocity  (ms-1) 

Well 

cell  voltage  (V) 

direction  along  the  channel  length 

y 

direction  normal  to  the  channel  length 

Greek  letters 

a 

excess  coefficient 

Q^area 

reaction  area  coefficient 

P02 

mole  fraction  of  oxygen  in  air  (20.9%) 

£h2 

mole  fraction  of  hydrogen 

Tf) 

overpotential  for  the  oxygen  reaction  (V) 

pi 

dynamic  viscosity  (N  s  m-2) 

P 

density  (kg  m-3) 

Pm, dry 

density  of  a  dry  membrane  (g  cm-3) 

0 

water  content  in  stream 

Gm 

membrane  conductivity  cm-1) 

Subscripts 

1A 

per  ampere 

air 

dry  air 

avg 

average 

A 

anode 

C 

cathode 

cell 

single  fuel  cell 

concentration  concentration  of  species  in  the  streams 

drag 

electro-osmotic  drag 

e 

electron 

H2 

hydrogen 

h2o 

produced  water 

in 

inlet  of  channel 

MW 

water  in  membrane 

n2 

nitrogen 

02 

oxygen 

OC 

open  circuit 

pressure 

partial  pressure  in  streams 

pump 

pump 

s 

stack 

sat 

saturation 

vapor 

water  vapor 

liquid 

liquid  water 

water 

all  water  including  vapor  and  liquid 

# 

cathode  or  anode 

cell  designs  with  coflow  and  counterflow  heat  exchangers. 
A  number  of  researchers  have  been  conducting  fuel  cell 
modeling  for  many  years  and  made  very  impressive  progress 
on  single  cell  modeling  [9-15].  These  models  emphasized 
important  characteristics  of  the  membrane  and  electrode 
as  well  as  a  detailed  description  of  the  water  content  in 
the  membrane.  To  the  authors’  knowledge,  most  of  those 
models  assume  that  there  is  no  temperature  and  pressure 
drop  along  the  channel.  Most  results  reported  were  at  high 
and  constant  stack  temperature  (e.g.  90-100  °C),  high  and 
constant  pressure  (e.g.  3  atm)  and  without  considering 
phase-change  effects.  But  in  practical  engineering,  the  PEM 
fuel  cell  is  usually  operated  at  varying  low  temperature  (e.g. 
65-75  °C)  and  varying  low  pressure  (e.g.  about  1.3  atm  for 
1  kW  portable  applications),  with  water  phase-change  inside 
the  fuel  cell  flow  channels. 

In  order  to  meet  these  challenges,  in  the  present  study, 
a  steady-state,  two-dimensional  mathematical  model  with 
pressure  effects,  water  phase-change  effects  and  detailed 
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(1)  Only  water  vapor  can  diffuse  into  electrode  and  pass 
through  the  membrane. 

(2)  The  electrode  layer  is  “ultra  thin”,  therefore  gas  diffusion 
through  the  porous  electrode  layer  is  neglected. 

(3)  The  gases  and  water  vapor  are  treated  as  being  well 
mixed,  the  mixture  is  assumed  to  be  an  ideal  gas. 

(4)  Water  is  produced  in  the  electrochemical  reaction  in 
vapor  form. 

(5)  Liquid  water  is  assumed  to  exist  in  the  form  of  small 
droplets  and  its  volume  is  negligible. 

(6)  No  voltage  drop  along  the  flow  channels  is  assumed. 

(7)  The  channels  in  the  single  cell  are  in  the  same  condition. 

(8)  The  solid  phase  temperatures  are  the  same  and  constant. 

2.2.  Mass  balance 


Fig.  1.  Schematic  diagram  of  PEM  fuel  cell  modeling  regions  (x-direction: 
perpendicular  to  the  paper). 

mass  and  heat  transfer  in  a  single  PEM  fuel  cell  was  devel¬ 
oped,  and  the  water  and  pressure  effects  on  a  single  PEM  fuel 
cell  performance  were  investigated  in  details. 

2.  Description  of  the  model 

A  typical  sandwich  construction  of  a  PEM  fuel  cell  is 
shown  in  Fig.  1.  The  model  regions  consist  of  a  proton 
exchange  membrane,  two  catalyst  layers  and  electrodes  at 
cathode  and  anode,  and  two  plates  with  flow  channels.  The 
v-axis  (normal  to  the  paper  in  Fig.  1)  is  parallel  to  the  gas 
channels.  The  temperature,  pressure  and  concentration  of 
gas  flow  will  be  predicted  along  this  direction.  The  y-axis 
is  normal  to  the  membrane.  The  hydrogen  ions  and  some 
water  molecules  transport  from  anode  to  cathode  along  this 
y-direction.  At  the  anode,  humidified  hydrogen  gas  is  sup¬ 
plied  at  a  given  excess  coefficient  into  the  flow  anode  chan¬ 
nels.  Hydrogen  gas  diffuses  through  the  porous  electrode  and 
reaches  the  anode  catalyst  layer,  then  the  chemical  reaction, 
2H2  — >  4H+  +  4e~ ,  releases  electrons  and  creates  hydrogen 
ions  which  pass  through  the  membrane  and  reach  the  cathode. 
At  the  cathode,  humidified  air  or  pure  oxygen  with  a  given 
excess  coefficient  flows  along  the  cathode  channels,  then 


The  amount  of  inlet  gases  is  calculated  according 
to  the  electro-chemical  reaction  for  a  PEM  fuel  cell, 
2H2  +  O2  =  2H2O  +  4e_.  One  equivalence  of  electrons  is 
1  mole  of  electrons  or  6.022  x  1023  electrons  (Avagadro’s 
number).  The  charge  on  one  mole  of  electrons  is  96,487 
coulombs  (C)  (Faraday’s  constant).  Thus,  the  charge  of  a 
single  electron  is  1.602  x  10-19  C.  One  ampere  of  current  is 
defined  as  1  C  s-1.  Therefore,  the  mole  number  of  electrons 
needed  per  second  for  1  A  of  current  is  as  follows: 

nQ  = - ^  =  1.03656546  x  10  5  mole  A  1  s  1 

6.022  x  1023 

(1) 

where  NEi  a  represents  the  number  of  electrons  for  1  A  and 
can  be  calculated  by  NEi  a  =  1/1.602  x  10~19.  Based  on  the 
reaction  equation,  the  theoretical  mole  numbers  of  consumed 
oxygen  and  hydrogen,  and  produced  water  for  1  A  current 
output  can  be  obtained  by  the  following  equations: 

no2,\A  =  \nQ  (2) 

^h2o,ia  =  (3) 

^h2,ia  =  (4) 

The  excess  coefficients  for  air  (oxygen)  or  hydrogen  is 

defined  as: 


actually  supplied  mole  number  of  air  (oxygen)  or  hydrogen 
theoretically  consumed  mole  number  of  air  (oxygen)  or  hydrogen 


diffuses  through  the  electrode  layer  and  meets  with  hydrogen 
ions,  the  reaction,  O2  +  4e_  +  4H+  — >►  2H2O,  occurs.  Water 
vapor  is  produced  along  the  cathode  channels  and  at  the  same 
time  electrical  energy  and  heat  are  released. 

2.7.  Basic  assumptions 


Therefore,  the  supplied  oxygen,  nitrogen  and  hydrogen  mole 
numbers  for  1  A  current  output  can  be  calculated  by  the  fol¬ 
lowing  equations: 


ft02,in,l  A  =  0'02^02,1  A 


«N2,in,l  A  =  «02,in,l  A  X 


1  ~  Po2 
Po2 


In  the  present  model,  the  following  assumptions  were 
employed: 


«H2,in,l  A 


0'H2^H2, 1  A 

^h2 
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where  po2  is  the  mole  fraction  of  oxygen  in  air  (/3o2  = 
20.9%)  and  /3h?  is  the  mole  fraction  of  hydrogen  in  anode 
(/3n2  =  1  in  the  present  study). 

For  generating  I  amperes  current,  the  molar  flow  rates  of 
oxygen,  nitrogen  and  hydrogen  for  the  single  channel  are 
therefore  evaluated  as  (mol  s-1): 


Nc,02,in  =  I  X 


^02,in,  1  A 
N ch 


due  to  electro-osmotic  drag  in  y-direction  is  evaluated  as 

[7]: 


N dragC}7)  — 


nd(x)I(x) 

F 


(15) 


where  I(x)  is  the  local  current  density  of  the  fuel  cell,  F  the 
Faraday’s  constant  and  nd  is  the  electro-osmotic  drag  coeffi¬ 
cient  representing  the  number  of  water  molecules  carried  by 
a  single  proton  calculated  by  [7] : 


nd(x)  = 


0.0049  +  2.02 a(x)  —  4.53 a2(x)  +  4.09 a3(x),  0  <  a(x)  <  1 

1.5849  +  0.159 (a(x)  -  1),  a(x)  >  1 


(16) 


Nc,N2,in  =  A^C,02,in  X 


1  ~  PO: 

Po2 


AA,H2,in  =  /  X 


^H2,in,l  A 
^ch 


(10) 


(11) 


The  components  of  mixture  vary  along  the  gas  channels 
and  the  local  molar  flow  rates  in  channel  are  evaluated  as 
follows: 


dNc,o2M  1{  .  j 

—  ^02,1A  X  I(x)  X  h  X  Of  area 

dx 

dNQN2(x)  _  0 

dx 


dNKn2{x) 

dx 


^H2,1A  X  I(x)  X  h  X  Of area 


(12) 

(13) 

(14) 


where  a(x)  is  the  activity  of  water  vapor  in  membrane.  Since 
the  membrane  is  placed  between  cathode  and  anode,  the  water 
vapor  activity  in  the  membrane  is  affected  by  the  water  vapor 
activity  in  the  flows  at  both  cathode  and  anode.  A  weighted 
average  water  activity  for  water  vapor  activity  in  membrane 
is  employed.  The  water  vapor  activity  at  anode  or  cathode, 
a#(x),  is  defined  similar  to  that  in  Reference  [7],  while  the 
pressure  p#(x)  is  a  variable,  not  a  constant  as  in  Reference 

[7]. 


a#(x)  = 


N#,  vapor  00 


X 


p#(x) 

P#,sat(x ) 


(/,  the  species  in  the  flow  stream  #) 


(17) 


The  diffusion  flux  caused  by  concentration  gradient  of  water 
can  be  written  as  follows: 


^concentration  (j)  —  ^MW 


(  9cmw 


(18) 


where  the  diffusion  coefficient  of  water  is  evaluated  as 

[7]: 


(0.0049  +  2.02a(x) 


Dmw  =  < 


1.59  +  0.159  [(a(x) 


4.53 a2(x)  +  4.09a3 (jc)  )D°  exp 


2416 


1  )]D°  exp 


2416 


1 


1 


303  Ts(x ) 


forO  <  a(x)  <  1 
for  a(x)  >  1 


(19) 


where  is  the  reaction  area  coefficient  that  accounts  for 
the  land  area  for  reaction  due  to  gas  diffusion  from  the  channel 
to  diffusion  layer. 

The  variation  of  vapor  and  liquid  water  along  the  channels 
is  more  complicated.  Vapor  water  transport  and  conden¬ 
sation  and  liquid  water  evaporation  are  considered  in  this 
model.  Three  water  transport  mechanisms  across  the  mem¬ 
brane,  according  to  Yi  and  Nguyen  [8],  are  considered  (a) 
electro-osmotic  drag:  since  hydrogen  ions  pass  through  the 
membrane,  the  water  molecules  are  carried  from  the  anode 
to  the  cathode;  (b)  back-diffusion  by  the  concentration  gra¬ 
dient  of  water:  because  the  water  concentration  is  different, 
some  water  molecules  diffuse  from  the  cathode  to  the  anode; 
(c)  convection  by  the  pressure  gradient:  water  moves  from 
higher-vapor-pressure  side  to  the  lower  one.  The  water  flux 


Convection  flux  caused  by  pressure  gradient  can  be  evalu¬ 
ated  in  a  similar  way  to  that  in  Reference  [8].  The  difference 
is  that  the  pressure  drop  along  the  channel  is  considered  in 
the  present  study  while  it  was  neglected  in  Reference  [8]. 


pressure 


cmw,cW  +  cmw,aW 
2 


kp 

x -  x 

MwW 


dp#,v  aporOO 

dy 


(20) 


where p#,y ap0r  is  the  water  vapor  pressure  in  the  anode  or  cath¬ 
ode  channels,  kp  the  permeability  of  water  in  the  membrane, 
pw(x)  the  water  viscosity  and  cmw,cW  and  cmw,aW  are  the 
concentrations  of  water  in  cathode  and  anode,  respectively, 
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and  are  expressed  as  [7]: 


/ 


cmw,#W  =  < 


\ 


Pm,  dry 

37  m,  dry 
Pm,  dry 

^Tm.dry 


[0.043  +  17.8a#(x)  —  39.8  a#(x)  +  36.0  a#(x)], 
[14  +  lA(a#(x)  —  1)], 


forO  <  a#(x )  <  1 
for  a#(x)  >  1 


(21) 


Pm, dry  and  Mm? dry  are  the  density  and  the  equivalent  weight 
of  a  dry  proton  exchange  membrane. 

Therefore,  the  change  of  water  flux  along  the  channel  at 
the  cathode  can  be  expressed  by: 


dAfc,  water  M 

dx 


—  [nu20  x  /(x)  +  2/drag(y)  ^concentration C)7) 

^pressure (j)]  x  hot  area  (22) 


where 

^7#,  vapor  00  —  ^7#,liquid(x) 

=  45070  -41.9[7#(x)-  273] 

+  3.44  x  10~3[7#(x)  -  273] 2 
+  2.54  x  10_6[T#(x)  -  273]3 
— 8.98  x  10~10[7#(jt)  -  273]4 


(26) 


The  water  flux  variation  along  the  channel  in  anode  can 
be  expressed  by: 


water  00 

dx 


[  3/ drag(}0  T  N concenatration  (y) 

T  ^pressure (j)]  Z  IlOi  area 


(23) 


The  mole  number  of  water  condensation  or  evaporation 
can  be  calculated  in  a  similar  way  to  that  in  Reference  [7], 
while  in  the  present  study,  the  pressure  is  different  at  different 
location  x  (pressure  was  constant  in  Reference  [7]): 


U#  represents  the  heat  transfer  coefficient  between  the  flow 
stream  #  and  stack. 

2.4.  Pressure  drop 

To  the  authors’  knowledge,  the  pressure  drop  of  mixture 
gas  in  the  fuel  cell  flow  channels  was  rarely  considered  in 
currently  available  fuel  cell  research  publications.  But  in 
industrial  design,  it  is  one  of  the  most  important  parameters 
simply  because  it  directly  affects  the  efficiency  of  a  fuel  cell 
system  and  is  directly  related  to  the  selection  of  the  system 
pump. 

The  saturation  pressure  (Pa)  can  be  expressed  in  terms  of 
the  local  temperature  [5]: 


/?sat(jc)  =  1.013  x  105  x  i0-2-1794+°-02953^W-273]-9-1837x10“5^W-173]2+1-4454x10^^W-273]3 


(27) 


dA^JiquidCO 

dx 


kchd  \  Af#,vaporOO 
RuT#(x)J  L £/##•/(*) 


X  P#(x)  -  p#,satC*0 


(/,  the  gaseous  species  in  stream  #) 


(24) 


Based  on  the  assumption  that  the  mixture  is  an  ideal  gas, 
local  volumetric  flow  (m3  s_1)  in  cathode  and  anode  can  be 
calculated  by  ideal-gas  law. 


Q#(x)  =  N#-'(x)  X  Ru  x 

i 


T#(x) 

P#(x) 


(/,  the  gaseous  species  in  stream  #) 


(28) 


2.3.  Energy  balance 

Based  on  the  assumption  that  the  temperature  of  solid 
phase  is  constant,  the  latent  heat  and  the  heat  transfer  to 
mixture  flows  from  solid  phase  are  considered  in  the  present 
study.  The  temperatures  of  the  mixture  flows  can  be  calcu¬ 
lated  as  [7], 


d7#(x) 


—  [Hw,  vapor  C*0 


#W,liquidM] 


dA^liquidC^) 

dx 


+  U#h(Ts  —  7#(x))  (/,  the  species  in  the  flow  stream#) 


Local  velocity  (ms  1 )  in  cathode  and  anode  can  be  cal¬ 
culated  as  follows: 

Q#(x) 

V#(x)  =  (29) 

where  A#  cross  is  the  cross-section  area  of  the  channel  in 
stream  #. 

Since  the  mole  fraction  of  gases  in  the  channel  varies,  local 
density  (kgm-3)  also  varies  due  to  different  components  in 
the  flow.  It  can  be  calculated  by: 


p#(x)  =  yf 

i 


N#j(x)  M#i 
£/  N#j(x)  x  1000 


P#(x ) 

T#(x)  x  Ru 


(25)  (i,  the  gaseous  species  in  stream  #) 


(30) 
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Local  dynamic  viscosity  can  be  calculated  by  interpola¬ 
tion.  /x/joo  is  the  gas  dynamic  viscosity  at  100  °C  and  /x;?o  is 
the  gas  dynamic  viscosity  at  0  °C.  The  temperature  scope  of 
flow  in  the  calculated  cases  is  0-100  °C,  so,  the  local  dynamic 
viscosity  is 


r#(x)  = 


X 


AM*) 

£,■  A r#,i(x) 

T#(X)  -  273 
100-0 


X  (Mz',1 00  -  Mi, o)  +  Ml, 0 


(/,  the  gaseous  species  in  stream  #) 


(31) 


where  7°  is  the  exchange  current  density  at  one  atmosphere 
of  oxygen  and  pc,o2  is  the  local  partial  pressure  of  oxy¬ 
gen  at  cathode.  It  is  affected  by  local  pressure,  which  is  not 
constant  in  the  present  study.  am(v)  is  the  membrane  conduc¬ 
tivity  and  calculated  by  the  following  equation  according  to 

[7]: 


am(x)  =  [0.00514  x  Mm,drycM :wO)  -  0.00326] 

Pm,  dry 


x  exp 


1268  x 


where 


(38) 


/ 


Cmw(*)  =  < 

< 


Pm,  dry 

A7rn,dry 
Pm,  dry 
A7m,dry 


[0.043  +  17.8a(v)  —  39.85 a2(x)  +  36.0  a3(x)], 
[14  +  1.4(a(x)  —  1)], 


forO  <  a(x)  <  1 
for  a(x)  >  1 


(39) 


The  average  current  density  is 


For  laminar  flow,  pressure  drop  in  each  control  volume 
can  be  expressed  as  (Pa): 


1  fL 

7avg  —  ~r  7  (v)  dx 

L  JO 


(40) 


dp#(x) 

dx 


=  P#(*)  x  f#(x) 


vj(x) 

2D 


(32) 


where  f#(x)  is  the  friction  factor  and  D  is  the  hydraulic  diam¬ 
eter  of  the  channel.  In  the  present  study,  the  channels  are 
assumed  to  be  straight,  thus  only  friction  loss  is  considered. 
The  local  pressure  is  calculated  by  the  pressure  at  inlet  sub¬ 
tracting  the  pressure  drop  from  the  inlet  to  current  control 
volume. 


p#(x)  =  p#t in  - 


d  p#(x) 

dx 


dx 


(33) 


Local  required  pumping  power  (W)  is: 


dP#  ,pump  (x)  =  Q#(x)  X  d p#(x) 


(34) 


Therefore,  total  required  pumping  power  (W)  in  stream  # 
used  by  the  designer  to  choose  a  proper  pump  to  maintain  the 
flow  can  be  given  by: 


pump 


=  Nch  X 


dp#(x) 

dx 


Q#(x) 


dx 


(35) 


3.  Solution  procedure 

The  model  equations  are  solved  by  numerical  method  at 
a  given  value  of  average  current  density,  7av g.  The  channel  is 
subdivided  into  n  control  volumes  of  equal  length  A L  =  Lin 
in  v-direction.  The  exit  values  at  the  Pth  control  volume  are 
the  inlet  values  at  the  (k+  l)th  control  volume.  A  set  of  dif¬ 
ferential  equations  is  replaced  by  algebraic  equations  based 
on  the  finite  difference  method.  Based  on  the  given  value  of 
7aVg,  the  flow  rates  for  hydrogen,  air,  water  vapor  and  liquid 
are  calculated.  Next,  the  guessed  value  for  the  cell  voltage, 
guessed  exit  temperatures  of  flow  streams  and  guessed  exit 
pressure  of  the  Pth  control  volume  are  chosen.  The  model 
equations  are  solved  to  get  a  set  of  current  densities  and  a 
new  /avg.  If  the  new  /avg  is  not  equal  to  the  given  7avg,  a 
new  cell  voltage  is  chosen  and  a  new  iteration  begins  until 
the  relative  error  between  the  calculated  7avg  and  the  given 
7avg  meet  the  convergence  criteria,  e.g.,  0.001.  Inside  the  iter¬ 
ations,  Newton-Raphson  method  [16]  is  used  to  solve  the 
non-linear  equations. 


2.5.  Cell  potential 


The  cell  potential  is  expressed  as  [7]: 


Lcell  —  Vqc  ??(*) 


I(x)tm 

crm(x) 


(36) 


where  Voc  is  the  open  circuit  potential  of  fuel  cell  and  r] (v)  is 
the  cell  activation  over  potential  and  is  calculated  as, 


RuTs(x)  ^  /  I(x)  \ 
0.5F  n  \  I  o  pc,o2(x)J 


(37) 


4.  Results  and  discussions 

For  convenience  of  discussion,  a  new  parameter  called 
relative  water  content  is  defined  as  follows: 

mole  number  of  water  (vapor  +  liquid) 

<fi#(x)  =  -  or 

mole  number  of  water  in  saturation 

,  ,  s  Af#  water (*)  /?#(*) 

mix)  =  — - - x - 

52i  N#AX )  P#,  satC*) 

(/,  the  gaseous  species  in  stream  #)  (41) 
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Table  1 


Size  of  channels  and  parameter  values  of  electrode  and  membrane 


Parameter 

Value 

Channel  number  of  cathode  and  anode 

4 

fiVch) 

Channel  length  (L) 

83.5  cm 

Channel  height  at  anode  ( d) 

0.1  cm 

Channel  width  at  anode  (h) 

0.1  cm 

Condensation  rate  constant  ( kc ) 

1.0s-1 

Dry  membrane  density  (pm,dry) 

2.0gcm~3 

Dry  membrane  equivalent  weight 

1 100  g  mol-1 

(Mm,dry) 

Membrane  thickness  (fm) 

0.01275  cm 

Diffusion  coefficient  of  water  in 

5.5  x  10-7  cm2  s-1 

membrane  ( D° ) 

Water  permeability 

1.58  x  10  14  cm 

Water  viscosity 

3.565  x  10-3  gcm-1s-1 

Fuel  cell  open-circuit  voltage  (Vqc) 

1.1  V 

Current  density  (/avg) 

1.0  A  cm-2 

Oxygen  exchange  current  density  ( I °) 

0.01  A  cm-2 

Heat-transfer  coefficient  (U) 

0.0025  W  cm-2  °C-1 

Calculations  were  performed  for  various  inlet  relative 
water  contents  on  anode  and  cathode  side,  respectively.  Pure 
oxygen  or  air  is  supplied  into  cathode  channels.  Channels 
are  in  rectangular  shape  and  the  size  of  the  channel  is  listed 
in  Table  1.  The  parameters  for  the  base  case  are  shown  in 
Table  2. 

Fig.  2  shows  the  effect  of  anode  inlet  water  content  on 
the  membrane  conductivity.  For  the  same  cathode  inlet  water 
content,  the  membrane  conductivity  increases  with  the  anode 
inlet  water  content  especially  near  the  inlet,  because  higher 
anode  inlet  water  leads  to  well-hydrated  membrane.  In  anode 
channel,  the  liquid  water  will  evaporate  to  generate  water 
vapor.  Fig.  3  shows  the  effect  of  cathode  inlet  water  con¬ 
tent  on  the  membrane  conductivity.  It  has  some  effect  on  the 
membrane  conductivity  near  the  inlet.  This  effect  decreases 
along  the  channel  because  the  vapor  pressure  increases  along 
the  cathode  channel,  and  condensation  occurs  when  the  vapor 
pressure  is  over  the  saturation  pressure. 

Fig.  4  shows  the  local  current  density  profiles  with  var¬ 
ious  anode  inlet  relative  water  contents  when  cathode  inlet 
water  content  is  fixed  at  1.0.  At  the  front  part  of  the  chan¬ 
nels,  the  membrane  is  well  hydrated,  which  leads  to  high 
membrane  conductivity  and  high  local  current  density.  The 
water  content  of  anode  side  dominates  the  membrane  perfor¬ 
mance.  Therefore,  the  current  density  increases  with  anode 

Table  2 

Parameter  values  for  the  base  case 


Parameter  Value 


Inlet  temperature  of  flow  at  cathode  and  anode 

333 

K 

Inlet  temperature  of  solid  phase 

333 

K 

Inlet  relative  water  content  at  anode 

0.5, 

0.75, 

1.0,  1.25 

Inlet  relative  water  content  at  cathode 

0.5, 

0.75, 

1.0,  1.25 

Inlet  pressure  at  anode  and  cathode 

1.1 

X  105 

Pa 

Excess  coefficient  of  flow  at  cathode 

2.0 

Excess  coefficient  of  flow  at  anode 

2.0 

Fig.  2.  Effect  of  anode  inlet  water  content  on  the  membrane  conductivity 
(pure  oxygen;  0c,in  =  1.0;  0a,in  =  0.5,  0.75,  1.0,  1.25). 


Fig.  3.  Effect  of  cathode  inlet  water  content  on  the  membrane  conductivity 
(pure  oxygen;  0a,in  =  1.0;  0c,in  =  0.5,  0.75,  1.0,  1.25). 


Fig.  4.  Effect  of  anode  inlet  water  content  on  the  current  density  (pure  oxy¬ 
gen;  0c,in=  1.0;  0a, in  =  0.5,  0.75,  1.0,  1.25). 
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Fig.  5.  Effect  of  anode  inlet  water  content  on  the  cell  voltage  (pure  oxygen). 

inlet  relative  water  content  near  the  inlet  up  to  the  inlet  with 
saturation  condition.  In  Fig.  4,  comparing  the  two  cases  of 
0a, in  =  1.0  and  0a, m  =  1 .25,  it  can  be  found  that  if  liquid  water 
is  supplied  (0a,in  =  1.25),  the  current  density  is  slightly  lower 
near  the  inlet.  This  is  because  the  change  of  membrane  con¬ 
ductivity  due  to  the  liquid  water  (0a,in  =  1 .25)  is  very  small 
near  the  inlet  (Fig.  2),  while  the  cell  voltage  for  0a,in  =  1.25 
is  higher  than  that  of  saturation  condition  (0a,in  =  1-0)  at  inlet 
(Fig.  5).  In  the  downstream  channel,  the  membrane  conduc¬ 
tivity  is  better  because  of  the  evaporation  of  liquid  water,  so 
the  current  density  is  slightly  higher  than  that  with  saturation 
condition  at  inlet.  If  the  anode  inlet  water  content  is  main¬ 
tained  at  1.0,  the  effect  of  cathode  inlet  water  on  the  local 
current  density  is  very  limited  (Fig.  6).  This  is  because  mem¬ 
brane  performance  and  water  transport  are  mainly  controlled 
by  the  water  content  of  anode  side. 

Fig.  7  shows  the  anode  inlet  water  effect  on  the  acti¬ 
vation  loss.  In  the  front  part  of  the  channel,  the  activation 


Fig.  6.  Effect  of  cathode  inlet  water  content  on  the  current  density  (pure 
oxygen;  0a,in  =  1.0;  0c,in  =  O.5,  0.75,  1.0,  1.25). 


Fig.  7.  Effect  of  anode  inlet  water  content  on  the  activation  loss  (pure  oxy¬ 
gen;  0c, in  =  1 .0;  0a, in  =  0.5,  0.75,  1 .0,  1 .25). 

loss  increases  with  anode  inlet  water  up  to  saturation,  which 
agrees  with  the  local  current  density  discussed  in  Fig.  4.  In 
the  rear  part  of  the  channel,  the  effect  is  reversed  due  to  the 
combined  effects  of  local  current  density  and  partial  pressure 
of  oxygen.  From  Fig.  8,  it  can  be  seen  that  the  activation  loss 
increases  with  the  increase  of  the  cathode  inlet  water  content. 
This  is  because  higher  inlet  water  in  cathode  leads  to  lower 
partial  pressure  of  oxygen.  The  cathode  inlet  liquid  water  has 
no  effect  on  the  activation  loss,  because  the  liquid  water  has 
no  chance  to  evaporate  in  this  case,  since  the  water  vapor 
in  the  cathode  is  always  saturated  or  over- saturated  due  to 
the  inlet  condition  with  liquid  water  and  water  production  in 
cathode. 

Figs.  9  and  10  show  the  effects  of  anode  and  cathode 
inlet  water  contents  on  the  ohmic  loss.  The  inlet  water  con¬ 
tents  of  both  sides  have  similar  effect  on  ohmic  loss,  i.e., 
increased  amount  of  inlet  water  leads  to  lower  ohmic  loss. 
This  is  because  the  increased  amount  of  inlet  water  leads  to 


Fig.  8.  Effect  of  cathode  inlet  water  content  on  the  activation  loss  (pure 
oxygen;  0a,in  =  1-0;  0c,in  =  O.5,  0.75,  1.0,  1.25). 
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Fig.  9.  Effect  of  anode  inlet  water  content  on  the  ohmic  loss  (pure  oxygen; 
0c,in  =  l.O;  0a, in  =  0.5,  0.75,  1.0,  1.25). 

a  higher  membrane  conductivity  as  shown  in  Figs.  2  and  3, 
which  dominates  the  ohmic  loss  although  the  higher  inlet 
water  contents  leads  to  higher  current  density  in  a  certain 
area  as  shown  in  Figs.  4  and  6. 

Fig.  1 1  shows  the  effect  of  anode  inlet  water  content  on 
the  net  water  flux  per  proton.  Near  the  inlet,  the  net  water 
flux  is  higher  due  to  high  electro-osmotic  drag  coefficient 
near  the  anode  inlet.  The  net  water  flux  increase  with  anode 
inlet  water  content  is  because  the  higher  anode  water  content 
provides  higher  electro-osmotic  drag  coefficient  as  shown  in 
Fig.  12.  The  cathode  inlet  water  content  has  almost  no  effect 
on  the  net  water  flux  in  this  case  as  shown  in  Fig.  13. 

Fig.  14  shows  the  effect  of  anode  inlet  water  content  on 
the  partial  pressure  of  hydrogen  and  oxygen  along  the  chan¬ 
nel.  The  higher  the  anode  inlet  water  is,  the  lower  the  partial 
pressures  are  for  both  hydrogen  and  oxygen.  For  hydrogen, 
higher  anode  inlet  water  leads  to  higher  partial  pressure  of 


Percentage  of  Length 


Fig.  11.  Effect  of  anode  inlet  water  content  on  net  water  flux  per  proton 
(pure  oxygen;  0C5in  =  1.0;  0a,in  =  0.5,  0.75,  1.0,  1.25). 
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Fig.  12.  Effect  of  anode  inlet  water  content  on  electro-osmotic  drag  coeffi¬ 
cient  (pure  oxygen;  0can  =  1.0;  0a,in  =  0-5,  0.75,  1.0,  1.25). 


Fig.  10.  Effect  of  cathode  inlet  water  content  on  the  ohmic  loss  (pure  oxygen; 
0a,in  =  LO;  0c,in  =  0.5,  0.75,  1.0,  1.25). 


Fig.  13.  Effect  of  cathode  inlet  water  content  on  net  water  flux  per  proton 
(pure  oxygen;  0a,in  =  1.0;  0c,in  =  0.5,  0.75,  1.0,  1.25). 
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Fig.  14.  Effect  of  anode  inlet  water  content  on  the  partial  pressure  of  hydro¬ 
gen  and  oxygen  (pure  oxygen;  0Cjin  =  1.0;  0a,in  =  0.5,  0.75,  1.0,  1.25). 

water  vapor,  which  decreases  the  partial  pressure  of  hydro¬ 
gen.  For  oxygen,  higher  anode  inlet  water  leads  to  higher 
water  vapor  transportation  across  the  membrane  that  leads 
to  lower  oxygen  partial  pressure.  Once  the  anode  inlet  has 
reached  saturation  condition,  its  effect  on  the  partial  pressure 
of  oxygen  is  very  small  because  the  water  vapor  of  anode  side 
is  enough  for  the  requirement  of  transportation  and  the  liquid 
water  evaporates  very  slowly  under  this  condition.  Therefore, 
the  increase  of  water  vapor  on  anode  side  because  of  liquid 
water  evaporation  almost  has  no  effect  on  the  pressure  of 
cathode  side.  Fig.  15  shows  the  cathode  inlet  water  content 
on  the  partial  pressure  of  oxygen  and  hydrogen.  Increased 
cathode  inlet  water  amount  leads  to  lower  partial  pressure 
of  oxygen  before  it  attains  saturation  condition  (0c,in=  1.0). 
When  (pc, in  >  1.0,  the  liquid  water  does  not  have  much  effect 
on  oxygen  since  the  liquid  water  in  cathode  channel  cannot 
evaporate  to  supplement  water  vapor  in  this  case;  because 
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Fig.  16.  Effect  of  anode  inlet  water  content  on  the  liquid  water  mole  number 
(pure  oxygen;  0c,in  =  1.0;  </>a,in  =  0.5,  0.75,  1.0,  1.25). 

the  partial  pressure  of  water  vapor  is  greater  than  the  satura¬ 
tion  pressure.  The  effect  of  cathode  inlet  water  on  hydrogen 
partial  pressure  is  very  limited. 

Fig.  16  shows  the  effect  of  anode  inlet  water  content  on 
the  liquid  water  mole  number.  It  is  clear  that  there  would  be 
no  liquid  water  in  anode  if  the  anode  inlet  water  content  is 
less  or  equal  to  one.  When  (p^-m  >  1.0,  in  anode  channel,  the 
liquid  water  decreases  along  the  channel  since  it  evaporates  to 
supplement  vapor  continuously.  In  cathode  channel,  the  liq¬ 
uid  water  increases  along  the  channel  since  vapor  condenses 
to  liquid  due  to  more  water  supplied  by  the  reaction.  The 
increased  anode  inlet  water  amount  corresponds  to  higher 
liquid  water  on  the  cathode  side  due  to  water  transport  across 
membrane  from  anode  side  to  cathode  side.  Once  the  anode 
inlet  reaches  saturation,  the  liquid  water  has  no  such  effect  as 
shown  in  Fig.  1 6  for  0a  in  =1.0  and  1 .25 .  As  shown  in  Fig.  17, 
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Fig.  15.  Effect  of  cathode  inlet  water  content  on  the  partial  pressure  of 
hydrogen  and  oxygen  (pure  oxygen;  0aan  =  1.0;  0c,in  =  0.5,  0.75,  1.0,  1.25). 


Fig.  17.  Effect  of  cathode  inlet  water  content  on  the  liquid  water  mole  num¬ 
ber  (pure  oxygen;  0ajin  =  1.0;  0c,in  =  0.5,  0.75,  1.0,  1.25). 
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Fig.  18.  Effect  of  anode  inlet  water  content  on  the  water  vapor  partial  pres¬ 
sure  (pure  oxygen;  0c,in  =  1.0;  0a,in  =  0.5,  0.75,  1.0,  1.25). 

higher  cathode  inlet  water  leads  to  higher  liquid  water  in  the 
cathode  channels  because  more  water  vapor  condensation 
occurs  for  the  cases  with  higher  cathode  water  content. 

Fig.  1 8  shows  the  effect  of  anode  inlet  water  content  on  the 
water  vapor  partial  pressure  along  the  channel.  Vapor  pressure 
decreases  along  the  anode  channel  due  to  water  transporta¬ 
tion  across  the  membrane  from  anode  to  cathode  side.  Vapor 
pressure  increases  along  the  cathode  channel  due  to  water 
transportation  and  electrochemical  reaction.  Increased  anode 
inlet  water  content  corresponds  to  increased  vapor  pressure 
for  both  sides.  As  shown  in  Fig.  19,  the  effect  of  cathode 
inlet  water  content  on  anode  vapor  pressure  is  relatively  small 
since  the  effect  of  back  diffusion  and  convection  due  to  water 
vapor  partial  pressure  difference  between  the  anode  and  cath¬ 
ode  is  correspondingly  small. 


Percentage  of  Length 

Fig.  20.  Effect  of  anode  inlet  water  content  on  the  Reynolds  number  at 
cathode  (air;  (f)c^n  =  1.0;  0a, in  =  0-5,  0.75,  1.0,  1.25). 

Reynolds  number  is  one  of  the  important  parameters  in 
fluid  flow.  Figs.  20  and  21  show  the  effects  of  anode  and 
cathode  inlet  water  content  on  the  Reynolds  number  at  the 
cathode  along  the  channel.  Reynolds  number  increases  in  the 
cathode  because  viscosity  decreases,  velocity  increases  and 
density  does  not  change  much  along  the  channel  as  shown  in 
the  authors’  previous  study  [17].  The  bigger  the  inlet  water 
content,  the  larger  the  velocity  and  Reynolds  number  at  the 
cathode.  For  the  present  study,  the  flow  is  obviously  in  lam¬ 
inar  flow  regime.  The  Reynolds  number  in  the  cathode  does 
not  change  when  0c?in  >  1.0,  this  is  because  the  liquid  water 
volume  is  neglected  in  the  present  calculations. 

The  pressure  loss  cannot  be  neglected  and  it  is  one 
of  the  important  parameters  for  fuel  cell  system  design. 
Figs.  22  and  23  show  the  local  pressure  in  the  cathode  along 
the  channel.  The  local  cathode  pressure  varies  almost  linearly. 
The  bigger  the  inlet  water  amount,  the  larger  the  pressure 


Fig.  19.  Effect  of  cathode  inlet  water  content  on  the  water  vapor  partial 
pressure  (pure  oxygen;  =  1.0;  0Cj in  =  0.5,  0.75,  1.0,  1.25). 


Fig.  21.  Effect  of  cathode  inlet  water  content  on  the  Reynolds  number  at 
cathode  (air;  0aan  =  1.0;  0Cjin  =  0.5,  0.75,  1.0,  1.25). 
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Fig.  22.  Effect  of  anode  inlet  water  content  on  the  local  pressure  at  cathode 
(air;  0c,in  =  l.O;  0a,in  =  O.5,  0.75,  1.0,  1.25). 

drop.  The  pressure  drop  in  the  cathode  is  around  7000  Pa, 
which  is  very  close  to  practical  operating  situation,  e.g.,  the 
Nexa™  Ballard  PEM  fuel  cell  operated  at  the  authors’  lab¬ 
oratory. 

Pressure  drop  for  different  control  volumes  is  not  equal. 
The  local  pressure  drop  increases  because  volume  flow  rate 
increases.  The  required  pumping  power  for  the  cathode 
increases  with  the  increase  of  channel  length,  which  is  shown 
in  Figs.  24  and  25.  The  greater  the  inlet  water  content  at  the 
anode  and  cathode,  the  larger  the  pumping  power  required. 
The  total  pumping  power  for  cathode  flow  of  a  single  cell  is 
proportional  to  the  number  of  channels  (four  channels  for  the 
present  study). 

Fig.  26  shows  the  effects  of  pressure  loss  in  the  flow  chan¬ 
nel  and  anode  inlet  water  content  on  the  power  of  a  single 
fuel  cell.  It  can  also  be  seen  that  when  the  pressure  loss  is 
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Fig.  24.  Effect  of  anode  inlet  water  content  on  the  required  pumping  power 
for  cathode  (air;  0C5in  =  1.0;  0a,in  =  0.5,  0.75,  1.0,  1.25). 


Fig.  25.  Effect  of  cathode  inlet  water  content  on  the  required  pumping  power 
for  cathode  (air;  =  1.0;  0c,in  =  0.5,  0.75,  1.0,  1.25). 


Fig.  23.  Effect  of  cathode  inlet  water  content  on  the  local  pressure  at  cathode  Fig.  26.  Effects  of  anode  inlet  water  content  and  pressure  loss  on  the  power 

(air;  0aan  =  1.0;  0Cjin  =  0.5,  0.75,  1.0,  1.25).  of  a  single  fuel  cell  (air). 
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Table  3 

The  effects  of  pressure  loss  and  anode  inlet  water  content  on  the  power  of 
the  single  cell 


0  a, in 

0.5 

0.75 

1.0 

1.25 

Voltage  (V) 

0.02 

0.13 

0.21 

0.25 

Current  (A) 

33.4 

33.4 

33.4 

33.4 

Power  (W) 

0.668 

4.342 

7.014 

8.350 

Power  loss  due  to  pressure  loss  (W) 

0.334 

0.334 

0.334 

0.334 

Percentage  of  power  loss  due  to 

50 

7.7 

4.8 

4 

pressure  loss 


considered,  the  cell  power  is  lower  than  that  without  con¬ 
sidering  the  pressure  drop.  The  most  important  reason  is 
that  the  pressure  loss  decreases  the  partial  pressure  of  oxy¬ 
gen,  which  leads  to  higher  activation  loss.  In  addition,  the 
pressure  loss  decreases  the  water  content  in  the  membrane, 
which  decreases  the  membrane  conductivity  and  increases 
the  ohmic  loss.  In  the  case  of  a  fixed  average  current  density, 
the  higher  the  anode  inlet  water  content  (0a,in  =  0.5-1.25), 
the  higher  the  fuel  cell  power,  i.e.,  the  fuel  cell  performance 
can  be  improved  by  humidifying  the  anode.  Further,  inject¬ 
ing  certain  amount  of  liquid  water  into  the  anode  inlet  could 
improve  the  fuel  cell  output  power.  Table  3  gives  detailed 
information  of  this  effect.  It  can  be  seen  that  the  percentage 
of  power  loss  due  to  the  pressure  drop  increases  with  the 
decrease  of  anode  inlet  water  content.  The  lower  the  fuel  cell 
power  is,  the  higher  this  percentage  would  be. 

5.  Conclusions 

The  model  for  reacting  fluids  on  both  the  air  and  hydrogen 
plates  of  a  single  PEM  fuel  cell,  with  water  transport  across 
membrane,  water  phase-change  effect  and  pressure  variation 
along  the  channel,  has  been  developed  and  the  predictions  of 
this  model  can  be  used  for  the  optimization  of  a  PEM  fuel  cell. 
In  addition,  the  effects  of  anode/cathode  inlet  relative  water 
content  on  fuel  cell  performance  were  investigated.  Based  on 
this  study,  the  following  conclusions  can  be  drawn: 

1 .  The  humidification  of  both  anode  and  cathode  sides  are 
very  important  for  the  performance  of  PEM  fuel  cells. 

2.  The  Water  content  of  the  anode  side,  which  tends  to  dry¬ 
out,  dominates  membrane  performance  and  plays  a  key 
role  in  the  performance  of  PEM  fuel  cells. 

3.  Proper  liquid  water  injection  at  the  anode  inlet  could 
improve  fuel  cell  performance  since  the  liquid  water  could 
evaporate  to  generate  water  vapor  and  thus  make  the  mem¬ 
brane  hydrated. 

4.  Liquid  water  injection  at  the  cathode  inlet  is  not  advo¬ 
cated  since  it  would  not  improve  the  cell  performance  but 
increases  the  burden  for  water  removal. 

5.  The  pressure  drop  in  the  PEM  fuel  cell  flow  channels 
increases  the  pumping  power  cost.  Therefore,  particular 
attention  must  be  paid  to  the  pressure  drop  when  design¬ 
ing  a  fuel  cell,  especially  if  air  is  used  rather  than  pure 


oxygen.  Pressure  loss  is  one  of  the  important  parameters 
that  affect  total  system  efficiency  and  optimization. 

6.  When  inlet  humidification  is  applied  to  improve  cell  per¬ 
formance,  the  increase  of  pumping  power  and  the  burden 
of  water  removal  because  of  humidification  should  also 
be  considered  in  the  system  design. 
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